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We have investigated nucleation of merons in double quantum dots when a lateral dis- 
tortion with a reflection symmetry is present in the confinement potential. We find that 
merons can nucleate both inside and at the edge of the dots. In addition to these merons, 
our results show that electron density modulations can be also present inside the dots. 
An edge meron appears to have approximately a half integer winding number. 
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1. Model Hamiltonian 

There is a range of values of the strength of the confinement potential for which 
spontaneous interlayer phase coherence is stable in double quantum dots^. Out- 
side this stability range the shape of the groundstate is expected to undergo a 
reconstruction^. It is the competition between the confinement potential and the 
Coulomb interaction which drives such a reconstruction. When the confinement po- 
tential weakens the repulsive Coulomb interaction pushes electrons from each other, 
and can create density depletions inside the droplet. 

The interplay between the reconstruction and the formation of merons in the 
presence of an external potential is unclear. We investigated within a Hartree-Fock 
(HF) theory how merons nucleate in the groundstate when a lateral distortion with 
a reflection symmetry is present in the confinement potential^ . We found that when 
a sufficiently large distortion is present merons can be nucleated. Both density 
depletions and merons can coexist in the reconstructed states. In the present work 
we demonstrate by integrating the electron density that these merons have total 
charge ±1/2. In addition we present numerical evidence for the presence of edge 
merons, which appear to have a half-integer winding number. 

An electron in each dot is in a parabolic potential ^m*n^r^. In the strong 
magnetic field limit the energy of an electron with the angular momentum m in 
a quantum dot is e„i — + hiUp{m +1), where lUc is the cyclotron frequency, 
hup = m*r2^^^, and £ is the magnetic length. Hereafter we will measure hojp in units 
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of e^/ei, using 7 — hujp/{e^ /e£) which is a dimensionless quantity proportional to 
the confinement potential strength (For GaAs 7 = 0.131{fujp[meV]f /{B[T]f/^). 
The real spins of all electrons arc assumed to be frozen aligning with the magnetic 
field. We consider a distortion in the confinement potential of the dot that may be 
generated by attaching a metallic stripe vertically on the cylindrical metallic gate 
surrounding the vertical quantum dots. The vertical distance between two dots is d. 
Let us make a simple model for this type of distortions. We take a separable form 
for the distortion potential U^^g{f,z) = vg{'r)h{z), where v is the strength of the 
distortion. The matrix elements of this potential are 

{niTi\Udis\n2T2) = {us6riT2+Ud{'^ - Srir2)) < ni\g{r)\n2 >, (1) 

where \njTj) = c\^.^. |0) denotes an electron with the angular momentum rij in the 
dot labeled by Tj. The parameters Ug and Ud are defined as 



dz/i(z) 



2 



2 

(2) 



Ud = v j dzh{z)cp* (^^ - V> (^^ + ■ (3) 

The parameter Us {ud) denotes the strength of the matrix elements of the disorder 

potential between intra (inter) layer electrons. The wavcfunctions (z + |) and 
^{z — I) are, respectively, the left and right dot wavefunctions along the vertical 
direction. 

The total Hamiltonian is H = Hq + V, where 

(4) 



ma ma m'a' 



and 

XI ^ {{mim2\Vs\m3m4) cl„^^cl„^„Cm3crCmiO- 
mimimsmi CT=Ti 

+ (mim2|Fd|m3m4> cj^^^cj^^ .^c^g -^Cm^,^^ . (5) 

The first term of Hq represents the confining energy due to the parabolic potential. 

In the second term, is the distortion potential. In the interaction Hamiltonian 
V , the term Vg {Vd) represents the Coulomb intca'action potential between the elec- 
trons in the same (different) dots. In the spatial coordinates, those interactions are 
given by Vs{r) = e^/er and Vd{r} = /s\f+dz\, respectively, where e is the dielec- 
tric constant. The spontaneous intcrlayer phase coherent state (maximum density 
droplet state"^'**), and meron excitations are well approximated by a HF theory^. 
The pseudospin density is computed with the calculated density matrix. 

2. Results 



Fig. 1 displays total electron densities for two different strengths of the lateral 
distortion with the center at the coordinate (—5, 5)£ by a Gaussian lateral potential 
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Fig. 1. Electron densities in units of l/H."^ for (us,u^) = (1.5, 0.25){e^/e£) and {us,u^) = 
(1.8,0.3)(e2/e£). 



with the range a — 3£. Results for 7 — 0.085 are shown for two different strengths 
of the distortion {us,Ud) — (1.5, 0.25), (1.8, 0.3)(e^/e£). The maximum density is 
bigger than l/27r^^, which imphes that the total local filling factor is bigger than 
one where the merons are located. This is consistent with the presence of merons 
with negative charges on the background of the total electron density l/27rf^. At 
{us,Ud) — (1.5, 0.25)(e^/e^) a meron is formed with the charge — e/2 and vorticity 
— . For the bigger value {us,Ud) — (1.8, 0.3)(e^/e^) two merons appear and the 
vorticities for the lower and upper merons are — and +, respectively. Integrating 
the electron density over the range that satisfy n(r) > ^^'-^ excess 

charges to be — 0.54e and — 0.98e, which indicates that one meron exists in the 
upper figure of Fig. 1 and a meron pair exists in the lower figure of Fig. 1. We 
have attempted to evaluate the total excess charge carried by the merons using the 
expression^ 

AiV = /d^r e^^m(r) • [9^m(r) x a^m(r)], (6) 

87r J 

where the polarization vector m(r) is f{r)/n{r) {f{r) — 2s /h and s is the pseu- 
dospin). However, this method gave inaccurate results. We believe this is because 
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the expression is valid only for slowly varying pseudospin fields. 




-5 5 

x/i 

Fig. 2. Arrows arc projections of r on the xy plane and the gray scale indicates the z-component 
of r. The dashed circle is the edge of the electron droplet and + and — indicate the core positions 
of merons with positive and negative winding numbers. The cross sign is where the impurity is 
located. The black dot represents the center of an edge meron, shown in the inset, and the small 
circles show electron density depletions. 

Fig. 2 displays how the pseudospin texture for 7 = 0.066 and {us,Ud) = 
(0.24, 0.04)(e^/rf). We observe merons inside the dots. The vahics of the polar- 
ization TO^(r) = Tz{r)/n{r) at the meron cores are, within numerical accuracy, +1 
and —1. This implies that the charges should be — e/2 since the vorticities of these 
two objects are opposite. The total charge of the merons is thus — e. Note that 
there are also five dc^nsity depletions inside the dots (the open circles in the fig- 
ure). These modulations of electron densities are not simply related to the shape of 
Upaiair) + Udis{r)- Since the densities are inhomogeneous the HF single particle po- 
tential must be included into consideration. It is non-trivial to anticipate the shape 
of the HF single particle potential since it reflects the delicate interplay between 
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Fig. 3. Magnified view of the spin texture of an edge meron with nearly a half-integer winding 
number. 



the Hartree self energy, exchange self energy, and confinement potential^. 

Wc also observe an edge meron near the edge in Fig. 2. A magnified view of the 
pseudospin texture is shown in Fig. 3. We draw a circle around the center of this 
object and see how the in-plane pseudospins rotate along the circle. The in-plane 
pseudospins nearly flips as the half circle is followed inside the dot. The magnitude 
of in-plane pseudospin values decreases to zero as the edge region is passed. We 
believe thus that the edge meron has approximately a half integer winding number. 
It is complicated to calculate accurately the excess charge of the edge meron since 
the electron density varies strongly near the edge even in the absence of an edge 
meron. It would be worthwhile to develop an analytical theory for edge merons. 
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